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Abstract 

We study in QCD the physics of deeply-virtual Compton scattering 
(DVCS) — the virtual Compton process in the large s and small t kinematic 
region. We show that DVCS can probe a new type of off-forward parton dis- 
tributions. We derive an Altarelli-Parisi type of evolution equations for these 
distributions. We also derive their sum rules in terms of nucleon form- factors 
of the twist-two quark and gluon operators. In particular, we find that the 
second sum rule is related to fractions of the nucleon spin carried separately 
by quarks and gluons. We estimate the cross section for DVCS and compare 
it with the accompanying Bethe-Heitler process at CEBAF and HERMES 
kinematics. 
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I. INTRODUCTION 



The Compton process, which refers to elastic scattering of a photon off a charged object, 
has played an important role in the history of Quantum Electrodynamics: It provided one 
of the early evidences that the electromagnetic wave is quantized, and hence has the nature 
of particles 0. The role of the Compton process in studying the structure of hadrons 
has been explored since 50's, when Low's low-energy theorems 0, analogous to the well- 
known Thomson cross section, are derived. Those theorems assert, for instance, that at 
sufficiently low energy the spin-dependent part of the Compton amplitude is determined by 
the anomalous magnetic moment of a composite system. Going to higher-order terms in the 
low-energy expansion, one finds the electric and magnetic polarizabilities 0. In recent years, 
experimental and theoretical works in measuring and understanding the polarizabilities of 
the nucleon and pion have flourished 0. 

Generally speaking, however, the Compton process on a composite system is quite com- 
plicated. When a point-like constituent absorbs an incoming photon, the system becomes 
excited. As it propagates in time, the system eventually emits a photon and comes back 
to the ground state. Quantum-mechanical propagation of a composite system is difficult to 
handle theoretically, except in spacial kinematic regions. The low-energy theorems exist for 
the Compton scattering on the nucleon because the intermediate propagation at low energy 
is dominated by the nucleon itself @. Another kinematic region known to have simple scat- 
tering mechanism is where the t-channel momentum transfer is large, i.e., the nucleon has a 
large recoil |J . In this case, perturbative Quantum Chromodynamics can be used to under- 
stand the intermediate propagation. In fact, according to the so-called the power counting 
rule ||, the most important intermediate states are those created from three valence quarks 
through hard-gluon exchanges. Simple as it may be, one still has to compute hundreds of 
Feynman diagrams to obtain the scattering cross section. 

The purpose of this paper is to study the Compton scattering by a virtual photon in 
a special kinematic limit. Assuming the virtual photon is generated by inelastic lepton 
scattering, we are interest in the Bjorken limit, i.e., the energy and momentum of the 
virtual photon going to infinity at the same rate. We shall call the process deeply-virtual 
Compton scattering (DVCS). As we shall discuss in the next section, the basic mechanism 
for DVCS is a quark absorbing the virtual photon, immediately radiating a real photon and 
falling back to the nucleon ground state. Thus the physics of DVCS is quite simple. 

Our interest in studying DVCS is generated from the fact that it offers a way to mea- 
sure the off-forward parton distributions (OFPDs), a new type of parton distributions that 
generalize the usual parton distributions and the nucleon form factors. When taking mo- 
ments of OFPDs, one gets form factors of the spin-ra, twist-two quark and gluon operators. 
When going to the forward limit, OFPDs become the usual quark and gluon distributions. 
Because the spin-two, twist-two operators are part of the energy-momentum tensor of QCD 
and because the form factors of the energy-momentum tensor contain information about 
the quark and gluon contributions to the nucleon spin, DVCS may provide a novel way to 
measure the fraction of the nucleon spin carried by the quark orbital angular momentum, a 
subject of great current interest @. 

The presentation of the paper goes like this. In Section II, we consider DVCS in QCD at 
the tree level. We identify the off-forward parton distributions from the Compton amplitude 
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and then study some simple aspects of the distributions, such as sum rules. In Section III, 
we study the leading-log evolution of the OFPDs. The results are presented in the form 
of generalized Altarelli-Parisi equations. In sections IV, we work out the cross sections for 
DVCS and the accompanying Bethe-Heitler process and their interference. To be complete, 
we consider all cases including unpolarized, double-spin and single-spin processes. Some 
estimates are given at the CEBAF and HERMES kinematics. The final section contains 
comments and discussions. 

A preliminary account of the DVCS process is discussed in a Letter paper by this author 
||. Subsequently, Radyushkin studied the scaling limit at A 2 = from a different angle ||. 
Actually, processes similar to DVCS were first considered by Geyer et. al. in studying the 



anomalous dimensions of light-ray operators |T(J , where the "interpolating functions" were 



introduced. The evolution of these functions is found to interpolate the Brodsky-Lepage 
and Altarelli-Parisi equations. Similar objects were also considered by Jain and Ralston in 
the context of studying the violation of helicity selection rule and the effects of transverse 
momentum [pTl . 



II. DEEPLY- VIRTUAL COMPTON SCATTERING AT LEADING ORDER 

This section is devoted to studying the virtual Compton scattering in deep-inelastic 
kinematics and at the leading order in perturbative QCD. The main result is that DVCS 
is dominated by single-quark scattering, and therefore the amplitude can be expressed in 
terms of the off-forward parton distributions. We also study sum rules of these distributions 
and show that the second sum rule is related to the total quark (gluon) contribution to the 
spin of the nucleon. 

We picture the virtual Compton scattering in Fig. 1, where a nucleon of momentum P^ 
absorbs a virtual photon of momentum q^, producing an outgoing real photon of momentum 
g'ti _ qH _ ^/x an( j a reco j} nucleon of momentum P' M = P M + A M . We focus on the deeply- 
virtual kinematic region of g M , namely, the Bjorken limit: Q 2 = —q 2 — > oo, P ■ q — > oo, and 
Q 2 / P ■ q finite. In this region, the quark that absorbs the virtual photon becomes highly 
virtual and hence propagates perturbatively. The simplest mechanism to form the Compton 
final state is for the quark to promptly radiate a real photon and fall back to the nucleon 
ground state. This "hand-bag" subprocess is shown in Fig. 2(a). 

In QCD, more complicated "tree" subprocesses are possible. By tree, we mean per- 
turbative diagrams in which every vertex is next to the nucleon blob. For instance, the 
highly- virtual quark can interact with gluon fields in the nucleon, as shown in Fig. 2(b); 
or it can transfer its virtuality to another quark through one-gluon exchanges, as shown 
in Fig. 2(c). A detailed calculation shows that both subprocesses are suppressed by l/Q 2 
relative to the hand-bag diagram, except when the polarization of the gluon in Fig. 2(b) 
is longitudinal. Fortunately, in the light-cone gauge A + = 0, the longitudinally-polarized 
gluons do not contribute by definition. Thus in the deeply- virtual limit, the single-quark 
process indeed dominates the Compton scattering. 

Of course, one can decorate the hand-bag with radiative loops, such as those shown 
in Figs. 2(d) and 2(e). Those diagrams in certain kinematic region produce leading-log 
corrections to the simple hand-bag, as we shall discuss in the next section. In other kinematic 
regions, they give rise to order-a s (Q 2 ) radiative corrections and hence can be ignored in 
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the deeply- virtual limit. One exception is the gluon vertex correction for the real photon 
in Fig. 2(e), where as the two quark lines have momenta predominantly parallel to the 
photon momentum, the diagram has an infrared divergence reflecting the non-perturbative 
photon wave function. Fortunately, a simple calculation shows that the contribution in 
this kinematic region is suppressed by l/Q 2 relative to the hand-bag due to the hard-gluon 
propagator. 

Thus, in the remainder of this section we concentrate solely on the dominant subprocess 
in Fig. 2(a). 

According to Feynman rules, the hand-bag diagram corresponds to the Compton ampli- 
tude, 

+^V (1 _ g) V, + ^M( t)} , (!) 

where /i and v are the polarization indices of the virtual and real photons and M{k) is a 
quark density matrix, 

M(k) = J e ikz d A z(P'\i){-az)i){{l - a)z)\P) , (2) 

where < a < 1 reflects the arbitrariness of the looping momentum k^. To proceed further, 
it is convenient to define a special system of coordinates. We choose and P M = (P+P')^/2 
to be collinear and in the z direction. Introduce two light-like vectors, p M = A(l, 0, 0, 1) and 
n M = (1, 0, 0, — 1)/(2A), with p 2 = n 2 = 0, p ■ n = 1, and A arbitrary. We expand other 
vectors in terms of p M , and transverse vectors, 

P"=p^+ (M 2 /2K , 

5" = + (Q 2 /2£K , 

= -f (p" - (M 2 /2K) + A^ , 
p = (k ■ n)p^ + (k ■ p)n" + kl , (3) 



where M 2 = M 2 - A 2 /4 and £ = (-P • q + y/(P ■ q) 2 + Q 2 M 2 )/M 2 . The £ variable is 
analogous to the Bjorken variable xb in deep-inelastic scattering and is bounded by and 
\J — A 2 / \J M 2 — A 2 /4. We neglect components of four vectors which do not produce large 
scalars in the Bjorken limit. Introducing the factor / dx^e lX<yX ~ k ' n ' > = 1 and integrating over 
k^ and we simplify the Compton amplitude to, 



T^(P, q, A)= ±(<T - p"n" - p v n») £ dx + 



x 



If (s, £, A 2 )f/(P') *tf(P) + E(x, £, ^)U{P') iaa "^ U{P) 



v a;-f/2 + ie x + £/2-ie / 



x 



A . T7 - 

#(*,£, A 2 )f/(P') fa 5 U(P) +E(x,£, A 2 ) w t/(P')75f/(P) 



(4) 
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where we have chosen a = 1/2 for symmetrical reason. H, H, E and E are off- forward, 
twist-two parton distributions defined through the following light-cone correlation functions, 

r d\ 

/ e iXx (P'M-\n/2)^(Xn/2)\P) = H(x, f, A 2 )U(P')YU(P) 



+E{x^A 2 )U{P')^^U{P) + ... 



— ( 

2vr 



-e 



< Xx (P'm-\n/2)Yl^(\n/2)\P) = H(x, £, A 2 )U (P')Yl 5 U (P) 



+E(x,^A 2 )U(P')^U(P) + ..., (5) 

where the ellipses denote higher-twist distributions. According to our definition, the initial 
nucleon and the active quark have the longitudinal momentum 1 + £/2 and x + £/2, respec- 
tively. [In covariant gauge, the longitudinal gluons produce a gauge link exp(—ig J X J n ■ 
A(an)da) between the two quark fields, restoring explicit gauge invariance of the light-cone 
correlations. Here we are working in the light-cone gauge n ■ A = 0, hence the longitudinal 
gluons and the gauge link never appear. Nonetheless, the result shall be taken as implicitly 
gauge-invariant . ] 

A systematic counting indicates that the virtual Compton scattering depends on twelve 
helicity amplitudes [|l^]. The above result implies that only four of them survive the Bjorken 



limit. An quick inspection shows all amplitudes with longitudinally-polarized virtual photons 
are sub-leading. For the transversely-polarized virtual photon scattering, the emitted photon 
retains the helicity of the incident photon at the leading order. This selection rule can serve 
as a useful check that deeply-virtual Compton scattering is indeed dominated by the single 
quark process. Of the four off-forward parton distributions, H and H conserve the nucleon 
helicity, while E and E flip the nucleon helicity. 

The off-forward parton distributions just defined have characters of both the ordinary 
parton distributions and nucleon form factors. In fact, in the limit of A M — > 0, we have 

If (x, 0,0) = q(x), #0,0,0) = Aq{x), (6) 

where q(x) and Aq(x) are quark and quark helicity distributions, defined through simi- 
lar light-cone correlations [0. On the other hand, forming the first moment of the new 



distributions, one gets the following sum rules Hill 



f^dxH{x,£,A 2 ) 


= F 1 (A 2 ) 


J l ^dxE(x,£,A 2 ) 


= F 2 {A 2 ) 


f^dxH(x,£,A 2 ) 


= G A {A 2 


j^dxE{x,i,A 2 ) 


= G P (A 2 



(7) 

where _Fi(A 2 ) and F 2 (A 2 ) are the Dirac and Pauli form factors and Ga{A 2 ) and Gp(A 2 ) 
are the axial-vector and pseudo-scalar form factors. [Usually the argument of form factors 
is the positive —A 2 . Here we omit the minus sign for simplicity.] 
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The second moment of the parton distributions is relevant to the spin structure of the 
nucleon. To see this, we first write down the angular momentum operator in QCD as the 
sum of quark and gluon contributions M, 



Jqcd — J q + J g i (8) 



where 



J q = I d 3 x x x T q 



d 3 x 

<3 



J g = J d 6 x x x (E x B) . (9) 

Here T q and E x B are the quark and gluon momentum densities, respectively. £ is the 
Dirac spin-matrix and D = d + igA is the covariant derivative. By an analogy with the 
magnetic moment, one can get the separate quark and gluon contributions to the nucleon 
spin if the form factors of the momentum density, or equivalently the energy-momentum 
tensor of QCD, are known at zero momentum transfer. Using Lorentz covariance and other 
symmetry principles, one can write down four form-factors separately for quark and gluon 
parts of the energy-momentum tensor, 

(P'\T^ g \P) = U(P')[A qtg (A 2 )^P v) + B qtg (A 2 )P^ia u)a A a /2M 

+C q , g {A 2 ){A»A u - g^A 2 )/M + C q , g {A 2 )g^M]u{P) , (10) 

where again P M = (P M +P M ')/2, A M = P^'—P^^ and U(P) is the nucleon spinor. Substituting 
the above into the nucleon matrix element of J qt9 , one finds fractions of the nucleon spin 
carried by quarks, J q , and gluons, J g , 

J q , 9 = \\KM + B q ^)} ; 

J q + J 9 = \- (11) 

According to the definition, the second moment of off-forward parton distributions yields 
the form factors of the energy-momentum tensor, 

f_ dxx[H{x, e, A 2 ) + E(x, £, A 2 )] = A(A 2 ) + B{A 2 ) , (12) 

where luckily the £ dependence, or C q (A 2 ) contamination, drops out. Extrapolating the sum 
rule to A 2 = 0, one gets J q ^ g . Note that only in this special application, we are interested in 
A 2 — > limit. In general discussions of DVCS, such limit is of course not necessary. 

By forming still higher moments, one obtains form factors of the twist-two operators of 
spin greater than 2. In general, there are many form-factors for each of the tensor operators; 
however, only special combinations of them appear in moments of the OFPDs. The relative 
weighting of the different form factors is determined by the variable £. 
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III. LEADING-LOG EVOLUTION OF OFF-FORWARD PARTON 

DISTRIBUTIONS 



In this section, we study the leading-log evolution of the off-forward parton distributions. 
Like the leading-log evolution of the usual parton distributions (Altarelli-Parisi equation 
||14||), there are many approaches to calculate it. Here we use the momentum-space Feynman 
diagram technique. For simplicity, our calculation is done in the light-cone gauge, although 
one is free to work entirely in covariant gauge. While the result of Subsection A has appeared 
in different forms in the literature before 0,|l^,|l^], the result of Subsection B is new. 

Evolution of the helicity-independent and helicity-dependent distributions are different, 
hence we treat the two cases separately. 



A. Evolution of Parton-Helicity-Independent Distributions 



In this subsection, we consider evolution of helicity-independent off-forward parton dis- 
tributions. We use a generic notation E s ^ NS (x, £, A 2 , Q 2 ) to denote singlet and non-singlet 
quark density, 



E 5 ^ 5 (x,e,A 2 ,Q 2 



Z7T 11 



(13) 



where the flavor indices and the gauge link have been ignored, and E G (x, £, A 2 , Q 2 ) to denote 
gluon distribution, 



E G (x,£,A 2 ,Q 2 



h I ^ Xx (n^ a (~n)F^n)\P)n,n u . 



(14) 



Since gluons are bosons, we have Ec{—x) = —Ec{x). The support for the light-cone variable 
x can be studied as in Ref. [F7| and is — 1 < x < 1. Since the evolution equations are 
independent of A 2 , we will omit the variable in the following equations. 

The evolution of the non-singlet quark density at x > £/2, where both quark lines 
represent quarks, takes the form, 



D Q E NS (x } ^Q 2 ) a s (Q 2 



where 



DlnQ 2 



d 



2tt 



1 —Pns(~, -)E NS (y, 6, Q 2 ) 

* y y y 



a s (Q 2 



-c f 



3 

2 + 



dy 



5/2 y — x — xe 



dy 



D\nQ 2 d\nQ 2 2tt 
The parton splitting function is calculated according to Fig. 3(a) 

x 2 + 1 - eii 



i/2 y — x — ie 



(15) 



(16) 



P NS (x,Z) = C F - 



l-x + ie){l-i 2 /A) 



(17) 



where Cp = 4/3 for SU(3) color group. The end-point singularity is cancelled by the 
divergent integrals in Dq/ D In Q 2 . Obviously, when £ = 0, the splitting function becomes 
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the usual Altarelli-Parisi evolution kernel. For — £/2 < x < £/2 where one of the two quark 
lines represents a quark and the other represents an antiquark, the evolution takes the form, 



D Q E NS (x,£,Q 2 ) _ a s {Q 2 



D\n Q 2 



2tt 



NS\~i I 



x y 



y y 



: dy P> ( x _^ 
-i y y y 



E NS (y,£,Q 2 ) , (18) 



where 



(19) 



When £ = 2, shifting the variable x — > x — 1 and then scaling the a; by a factor of 2, one finds 
that the evolution equation becomes the Brodsky-Lepage evolution equation for the pion 
wavefunction [18[]. For x < — £/2 where both quark lines represent antiquarks, the evolution 
takes the same form as Eq. (|15|), apart from the replacement f x — > — J^ 1 . 

The evolution of the singlet-quark density mixes with that of the gluon density. For 
x > £/2, the coupled evolution takes the form, 



D Q E s (x,£,Q 2 ) 

DlnQ 2 
D G E G (x,Z,Q 2 ) 

D\nQ 2 



a s (Q 2 )_ r^dy 
Jx y 



2vr 

a s (Q 2 ) 
2tt 



Pss(-, -)E s (y, £, Q 2 ) + 2n F P SG (-, i)E a {y, £, Q 2 ) 

y y y y 



• y 

+ Pgg 



pgsC-i -) x ~ (My, e, q 2 ) - E s (-y, e, q 2 : 

y y 2 v 

X £ 



y y 



Eg(v,€,Q 



where rip is the number of quark flavors and 



D 



G 



d 



a s {Q 2 

D\nQ 2 dlnQ 2 2vr 

The evolution kernels are, 



-a 



ii 



Up 



6 3Cj 



+ 



dy 



+ 



i/2 y — x — xe 



dy 



-f/2 y — x — xe 



(20) 



(21) 



x 2 + (1 - x) 2 - £ 2 /4 



P G5 (x,0 = Cf 



(1 - £ 2 /4)^ 



1 + fl 



£ 2 /4 



x(i-e/4) 

(x 2 - £ 2 /4) 



x(l-e 2 /4) ; 



2(l-x)(l+x 2 ) l+x-£ 2 /2 



r 



e 2 /4 



1 — x + ie 



(22) 



where Tp = 1/2 and = 3. As £ — ► 0, we again get the usual Altarelli-Parisi evolution 
kernels. 

For — £/2 < x < £/2, we have 



D Q E s (x,e,Q 2 ) _ a s (Q' 



D\n Q 2 



2tt 



1 dy x £ 

* J ss{ ?y 



^ dy . . x £ , 
12/ 2/2/ 



^5(2/, e,Q 2 ) 
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D G E G (x,C,Q 2 



» 2/ y y 



-p. 



i y 



SG\~ ) _ , 
1/ 1/ 



£ G (y,£,Q 2 



2tt 



1 / X ^ 



« S (Q 2 ) 

2vr 



x y 

x 
i 



y y 



1 dy 0/ ,z 



-i y 



y y 



x y 



y y 



, X 



P G G\~i > I Pnni—, — 



-1 y 



GG\~i 



y y 



where the evolution kernels are, 



x + £/2)(l-2x + £/2) 



P' SG (x,0=T F - 
P GS (x,0=C F 

P'ggM-- 



C- 



e(i + e/2)(i-e 2 /4) 
> + e/2)(2-x + e/2) 
<(i + e/2) 

(x 2 - e 2 /4) 

<(i-e/4) 



E G (y,£,Q 2 ), (23) 



2(1 + x 2 ) 



1-z + ie (l + f/2)(x-£/2) 



(24) 



For x < — C/2, the evolution equations are the same as Eq. (p0|), except /J 



B. Evolution of Parton-Helicity Dependent Distributions 



The helicity-dependent distributions are defined as follows: The singlet and non-singlet 
quark distributions are 



E StNS (x,£,A 2 ,Q 2 
The gluon distribution is, 

E G (x,^A 2 ,Q 2 ) = - 



^ Xx {Pm~n) M{\n)\P) . 



-_L f — e ^(P'\F^(--n)F v J-n)\P)n u n v , 
2xJ2n x 1 2 2 / m ^ ' 



(25) 



(26) 



where P a/3 = \e a ^ 6 F 1& . It is easy to see that E G (-x) = E G (x). 

The evolution of the non-singlet helicity-dependent quark density is exactly the same as 
that of the non-singlet helicity-independent quark density. 

For the singlet evolution, we consider mixing between the singlet quark and gluon dis- 
tributions. For x > £/2, the evolution equations are, 



D Q E s (x,^Q 2 ) = a s (Q 2 ) f dy 
I) In Q 2 2- Jx y 



AP SS (-, i)E 8 {y, e, Q 2 ) + 2n F AP SG (-, —)E G (y, £, Q 2 ) 

y y y y 



D G E G (x,^Q 2 ) a s (Q 2 ) f dy [ A ^ £ 1 , ~ 2 ~ , n 

n1 n2 — = — ^ — / — A.Ffcs(- - x o [ E s{y,€,Q ) + Es{-y,Z,Q 

DmQ 2 2ti Jx y y y 2 v 



+ AP GG (-,i)p G ( 2/ ,e,g 2 ) 

y y 



(27) 
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where the splitting functions are, 



AP ss (x,0=Pns(x,0 , 

x 2 -(l- x) 2 - e 2 /4 



AP SG (x,0=T F - 
AP GS (x,0 = C F 
AP GG (x,£)=C A 



;i-e/4) 2 



i-fi 



£ 2 /4 



x(i-e/4) 
(x 2 - e/±) 



xii-ei^y 



Ax(l-x) l + x-£ 2 /2 

1 + —^7 + 



x 2 - £ 2 /4 



1 — x + ie 



(28) 



Again, when £ = 0, the splitting functions are the usual spin-dependent Altarelli-Parisi 
kernel. For — £/2 < x < £/2, the evolution equations are, 



D Q E s (x^ } Q 2 ) a s (Q 2 



Pin Q 2 



D G E G (x,£,Q 2 
Pin Q 2 



2tt 

« S (Q 2 



dy 



,x i 



-^APU-, -) = - / — AP' SS (-, -1 



* 1/ 



2tt 
2tt 



2tt 



-2n, 



SSK ?y 
* y SGK y y' 



,x £ 



dy 



,x 



i y 



y y 



E s (y,Z,Q 2 



" ^AP' (-,-- 

i 2/ SG V y' 



E G (y,Z,Q 2 



x y 

x 



Gsv yy 



— AP'nJ-, -) - I —AP'nJ-, -- 



dy 



,x 



|(£s(y,£,g 2 )+£sH/,f,Q s 



1 —AP GG (-, -) 

* y GG y y 



-i 2/ GG Y 2/ 



P G (?/,£,Q 2 ) , (29) 



where the splitting functions are, 

*P' ss (x,£) = P' NS (x,t) , 

(x + £/2)(-l + £/2) 



AP^(x,£)=7> 
AP G5 (x,£) = CV 



e(i+£/2)(i-e/4) 
<(i + e/2) ' 



e 2 /4) 



4;r 



x + «e (l + £/2)(z-£/2) 



(30) 



For x < — £/2, the evolution equations are the same as Eq. (|2"?D, except f% — > — J^. 

Of course, in a reasonable range of Q 2 , the OFPDs do not evolve dramatically. Thus, a 
second check on the DVCS dynamical mechanism is to find small Q 2 dependences of relevant 
scaling functions. 



IV. CROSS SECTIONS AND ESTIMATES 

In this section, we calculate the cross section for electroproduction of a real photon off a 
nucleon. As shown in Fig. 4, we use k = (u, k) and k' = (a/, k') to denote the four-momenta 
of the intial and final electron, P = (M, 0) and P' = (E', P') the intial and final momenta 
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of the nucleoli, q' = (v',q') the momentum of the final photon. The differential cross section 
in the laboratory frame is, 

where M is the nucleon mass and T is the T-matrix of the scattering. Integrating over the 
photon momentum, we find, 

1 .^..nUj'duj'dQ e / „ , . .. 9 , d 3 P' , . 

= i T| Sp# M «* + F " * " F » »2iw ' (32) 

where (^-function reflecting the photon on-shell condition, which constrains the direction and 
magnitude of the momentum of the recoiling nucleon, P', 

s + M 2 - 2 {[y + M)E' -q-P')=0, (33) 

where s = (q + P) 2 , q^ = {y,q) = — Thus the phase space of the recoiling proton 
is specified by the solid angle d£lp>. Note, however, that for large s and v, there are two 
solutions of \P'\ corresponding to one orientation of P'. Physically one of them represents 
the recoil proton at the backward angle in the center-of-mass frame. Integrating over the 
magnitude of P', we get, 

1 P' 2 

d(J = — ^^u'du'dQe'dnp,-— — — -\T\ 2 , (34) 

32(2tt) 5 cjM \P'(u + M) -qE> cos <f)\ l 1 ' K J 

— * — * 

where is the angle between q and P', and the sum over two possible solutions of \P'\ is 
implicit. We choose the z-axis to be the direction of the incident electron and the x-axis in 
the plane formed by the initial and final electron momenta. In this coordinate system, the 
final electron has the polar angle 9. The polar and azimuthal angles of P' are denoted by 
Opt and (ppi, respectively. Note that \T\ 2 has the dimension of a cross section. 
Alternatively, we can integrate out the momentum of the recoil proton, 

1 v' 

da = — ; — r-= u' du 1 dVL e tdVt a i- :\T\ 2 , (35) 

32(27r) 5 c;M q \v + M - q cos <j/\ 1 1 ' K J 

where 0' is the angle between q and q'. The polar and azimuthal angles of q' are denoted by 
9 q i and <p q >, respectively. The constraint between the energy and direction of the outgoing 
photon is, 

s-M 2 -2u'(u + M -qcos(f)') = . (36) 



A. T-matrix 

We calculate the Feynman diagrams shown in Fig. 4. We assume for the moment that the 
scattering lepton is negatively charged (e > 0). The T-matrix for the Compton scattering 
part is, 
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7i = - e ^u(k')Yu(k)-T^* , 



(37) 



where u, u are the spinors of the lepton and e is the polarization of the emitting photon. 
is the standard Compton amplitude, 



T, u = iJ e-^ z (P'\TJ^z)J u {U)\P)d A z 



(38) 



In the deeply-virtual region, T^ v is expressed in terms of off-forward parton distributions in 
Eq. (4). 

The T-matrix for the Bethe-Heitler (BH) process is, 



T 2 = -e 3 u(k') 



— 7 M + Y 



-f 



u{k) — {P'\J^)\P) , (39) 



ft- A - ™ e + it ' ' ' ft'+ A - m e + it 

where m e is the mass of electron and will be ignored for the following discussion. The elastic 
nucleon matrix element is, 



(PU(0)|P) = U(P') k^(A 2 ) +F 2 (A 2 )^^I U(P) , 



(40) 



where U, U are the nucleon spinors and F\ and F 2 are the usual Dirac and Pauli form factors 
of the nucleon. 

The total T-matrix is the sum of the two above, 



T = T 1 + T 2 . 



(41) 



B. Unpolarized Scattering 

First, we consider the scattering process without polarizations. The square of the T- 
matrix has three terms. First, the pure Compton process gives, 



I^r = -Vvcww, 



(42) 



where the lepton tensor is, 



e$> = 2(W + k u k'^ - g» v k ■ k') . 



The hadron tensor is calculated from Eq. (4), 



vv vc 



J dxa(x) J dx'a*{x') 



(43) 



[-e(H(x, A 2 ) + E(x, £, A 2 ))(H(x', A 2 ) + E(x', A 2 )) 
+4H(x, C A 2 )H(x>, C A 2 ) - ^E(x, £, A 2 )E(x', £, A 2 )] 

+ J\x(3(x) J dx'(3*(x')[-e(H(x,^A 2 )+E(x,^A 2 ))(H(x\^A 2 ) + 



+4H(x, £, A 2 )H(x', £, A 2 ) + eE(x, £, A 2 )£(a;', £, A 2 ) (1 



4M 2 



)] 



(44) 
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where a(x) and (3(x) are, 



a{ - %) x-i/2 + ie + x + i/2-ie ' 

P(x) = : ^- . (45) 

Thus, for H(x) and E(x) less singular than x~ 2 and H(x) and E{x) less singular than 
the integrals are convergent. To simplify the expression, we have left out as usual the sum 
over quark flavors weighted by the electric charge squared. To find the product of the two 
tensors, one shall express the null-vectors n and p in terms of physical vectors P and q by 
inverting Eq. (3). The result is that |7^| 2 goes like l/Q 2 at large Q 2 , like the inclusive 
deep-inelastic process. 

The pure Bethe-Heitler process gives, 

\%\ 2 = J^WkHG-o , (46) 

where the lepton tensor is 

tfgg = 8 A , 4 [k ■ (k' + A)](2W + AW + AW - g^k' • A) 

+ - —r-: [k' ■ (k - A)](2Fr - AW - A u k^ + g^k • A) 
(k — A) 4 

\[2(k' + A) • (k - A) + A 2 \{k'W + h!W) + 2{k'WA ■ k 



(k - A) 2 (k' + A) 2 - 

kWA ■ k') + [(Jfe" - fc'^A" - (A; ,iy - fc^A") + (A;' - fc) • A^]A; • fc'l . (47) 



As A M — > 0, it diverges quadratically as expected from low-energy theorems. The hadron 
tensor is, 

= (Fx (A 2 ) + F 2 (A 2 )) 2 (A V - A" A") 

+(P + ^Y(P + f r4(F x (A 2 ) 2 - t^^(A 2 ) 2 ) , (48) 

which is well-known from elastic scattering. Although the product of the two tensors is 
complicated in general, it simplifies to 

21 A» [k-A k'-A) K ' 

as A 11 — > 0. Thus, for small A M , the cross section of the Bethe-Heitler process dominates 
that of DVCS. To have a clear DVCS signal, one must have a A 2 reasonably large (at least 
on the order of the nucleon mass) and a Q 2 not too large. 

To appreciate the relative and absolute sizes of the DVCS and BH cross sections, we 
shown in Fig. 5 some calculations at the electron beam energies u = 6 GeV (CEBAF after 
upgrading) and u = 30 GeV (DESY). For the Bethe-Heitler process, the cross section can 
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be evaluated accurately using the experimentally-measured nucleon form factors. For the 
DVCS cross section, we need a model for the OFPDs. Since at the moment we are interested 
in only a rough estimate, we assume E = E = and 



H(x,£,,A 2 ) = H(x,£,A 2 ) = q(x)e A2 ^ . 



(50) 



In Fig. 5a, we have shown the DVCS (solid) and BH (dashed) cross sections at incident 
electron energy uo = 6 GeV and scattering angle 9 = 18°, with a virtual photon v — 3 GeV 
and Q 2 = 1.76 GeV 2 . The recoil proton is detected in the electron scattering plane and at 
the same side of the scattered electron. In Fig. 5b, we have shown similar cross sections at 
uo = 30 GeV and scattering angle 6 = 3.5°, with the photon energy v — 7 GeV, Q 2 = 2.6 
GeV 2 . Since the BH cross section is calculable, experimental cross sections tell us about the 
DVCS cross section plus the interference cross section which we now turn to. 
The interference contribution has the following structure, 



T*T 2 + 7iT 2 * = -2^- 2 (&^ a ReH (pu)a + M a ReH [Ha ) , 



(51) 



where the lepton tensor has been divided into symmetric and antisymmetric parts according 
to the indices ji and v. The symmetric lepton tensor is 



(k-Ay 



(k ,fi k u + - A) a + [k'^k - Ay + k w {k - Ay\k a 

+ {k'»g ua + k' v g» a ){k- A) 

-g^[k ■ k'(k - A) a + k'-(k- A)k a + k ,a (k ■ A)] 

(k^k u + k w wyk' + Ay + [k^(k' + Ay + k v (k' + Ay\k ,c 

— {k^g va + k u g fia )(k' ■ A) 

-g» v [k ■ k'(k' + A) a + k ■ (k' + A)k' a - k a (k' • A)]] , 



(52) 



and the antisymmetric lepton tensor is 

2 

(k - A) 2 

-k' a (k»A v - k u A») 



(k ■ k')(A»g va - AV°) + {k' ■ A){k^g ua - k v g^) 
2 



+ 



(fc-jt')(Ay a -Ay a ) 



{k! + A) 2 

— {k ■ A)(k'^g va - k' v g^) + k a {k'»A v - k' v A»)] . (53) 

As A^ — > 0, the symmetric part dominates over the antisymmetric part. 

The hadron tensor is also separated into two contributions. The symmetric part is, 



H (^)a = _ p „ n u _ p u n ^ J 1 dxa *^ 



(Fi(A 2 ) + F 2 (A 2 ))(E(x, £, A 2 ) + H(x, £, A 2 )) [|a° + ^-n a ] 



+2P a [H(x,£,A 2 )F 1 (A 2 ) 



,2n 

A 2 
AM 2 



£(*,£, A 2 )F 2 (A 2 )] 



(54) 
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and the antisymmetric part is 

# [H° = - € ^n pP<T e^ 5a n 01 A s £ ^*(x)(F 1 (A 2 ) + F 2 (A 2 ))tf (:r, £, A 2 ) . (55) 

The final expression for products of the tensors above is quite lengthy and is omitted. It is 
generally believed that the interference contribution has a size between the BH and DVCS 
cross sections. Thus, if the BH cross section is much larger than that of DVCS, the DVCS 
amplitude might be accessible through the interference term which, for instance, can be 
extracted by comparing electron and positron scatterings, or by a direct subtraction. 



C. Double Spin Process 

In this subsection, we consider scattering with both lepton beam and nucleon target 
polarized. The structure of the cross section is similar to that of the unpolarized case. 
Indeed, the spin-dependent part of the T-matrix squared has three contributions. The pure 
virtual Compton process gives, 

\Ti\ 2 = -p^WyciH ■ (56) 

The spin-dependent lepton tensor is antisymmetric, 

= -2\ie^k a k' , (57) 

where A = ±1 represent the positive or negative helicity of the scattering lepton. The 
spin-dependent hadron tensor is also antisymmetric, 

Wfc = JdxJ dx'Re{a(x)P*(x')][4H(x,Z,A 2 )H(x',Z,A 2 )(S-n)(l - |) 

-2E(x, e, A 2 )H(x', A 2 )[(S ■ n)(£ - ^) - ^^(1 + |)] 

-(H(x, £, A 2 ) + E(x, £, A 2 ))E(x', £, A 2 )X [( s . A )(l + |) + (S ■ n)^] 



+E(x,^A 2 )E(x',^A 2 )^(S-A)\ , (58) 

where is the polarization of the nucleon with normalization S 2 = —M 2 . Using 
el ivap e is _ _2(g<xigPS — g a5 g^)^ one can straightforwardly work out the product of the 
two tensors. 

The pure Bethe-Heitler process gives, 



W = p l ^W mH . (59) 



The antisymmetric lepton tensor is 

£M = -8Xie^ al3 A a 

x - k' p {k' + A)-k + k p {k-A)-k' _ {k + k') p {k-k') j Q 



(A;' + A) 4 (k - A) 4 (k - A) 2 {k' + Af 
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and the antisymmetric hadron tensor is, 



ie^A a \2(F 1 (A 2 ) + P 2 (A 2 )) (F X (A 2 ) + ^^(A 2 )) S i 



+P 2 (A 2 )(P 1 (A 2 )+P 2 (A 2 ))i^P, 



(61) 



As A M — > 0, the lepton tensor is sub leading relative to the spin-independent counterpart 
shown in Eq. (47). Thus the spin asymmetry in the BH process vanishes in such a limit. 
Finally we consider the interference contribution, 

where the spin-dependent lepton tensor has both symmetry and antisymmetric parts. The 
symmetric part is, 

2\i 



T*T 2 + = -2^3 (A£^ a Re[AH { , u)a } + M^ a Re[AH [Ha ]) , 



(62) 



A& v)a = - \,J k' p e uapa + k! v t mpa - g^e Xapa k' x )A p k a 
(k — A) 2 v ' 

f ,, 2X \^ (k p t va(XJ + k v e papa - g pu e Xapa k x )A p k' a , 



(k 1 + A) 2 
and the antisymmetric part is 
2A 



(63) 



MiHa = (k-A\^ VfXJ ( {k ~ Arkp + ka{k ~ A)p + ^ ' A ) k '° 



At 



(64) 



(k - A) 2 

- (fc /+ A A) 2 ^ 0' + A ) QA; P + + " 

As A^ 0, the antisymmetric part dominates over the symmetric part. 

The spin-dependent hadron tensor is separated into two parts accordingly. The symmet 
ric part is, 

AH^ a = -i(g pu - p p n v - p u n p ) J* dxa*(x) 

"(F X (A 2 ) + F 2 (A 2 ))(E(x, £, A 2 ) + H(x, £, A 2 )) e" xp °n x S p A a 



+F 2 (A 2 )(H(x,£ t , A 2 ) + E(x,£, A 2 )) El e x ^p x S p n a A T 



-(F 1 (A 2 ) + F 2 (A 2 )) J E( a ;,e,A 2 )— e^P A ^A CT 
and the antisymmetric part is, 

Ai/H Q = i^^npPa J 1 dx(3*(x) 

(Pi (A 2 ) + F 2 {A 2 ))H{x, £, A 2 ) [(5 • n)2P a - £S a - (5 • A)n a 
-2PT(:r,£,A 2 )P 2 (A 2 )P a 



(65) 



-(P X (A 2 ) + P 2 (A 2 ))P(x, £, A 2 )^[(^ • A)P° + ^1 



2M 2 



P 2 (A 2 )P(x,£,A 2 )(S-A)P Q 



(66) 
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The product of tensors leads to a lengthy expression which we again omit. 



D. Single-Spin Process 

The virtual compton amplitude is complex because of the intermediate quark propa- 
gation. This gives rise to the so-called single-spin asymmetry, that is, the cross-section 
asymmetry depending on a single polarization. For instance, suppose the lepton beam is 
polarized and the hadron is unpolarized, the single-spin asymmetry is proportional to 

A e L ~ -2^ (A&^ a Im# (HQ + A£^ a hnH^ ]a ) . (67) 

On the other hand, if the nucleon is polarized and the lepton is unpolarized, the single-spin 
asymmetry is proportional to 

< T ~ -2^- 2 (£^ a Im[Atf (Ha ] + fH« l m [AH [Ma }) . (68) 

The size of these spin asymmetries directly reflect the relative contributions of the DVCS 
and BH processes. Their measurement and interpretation are very interesting although they 
only depend on the parton distributions at x — ±£/2. 



V. COMMENTS AND CONCLUSIONS 

In this paper, we studied some basic aspects of deeply-virtual Compton scattering. The 
motivation is that in the deeply-virtual kinematics, the scattering mechanism appears to be 
simple, and hence, one can learn some structural information from the process. The QCD 
analysis shows that it is the off-forward parton distributions that are probed. In a particular 
experiment, of course, one has to decide whether one is at the deeply- virtual kinematics. As 
we discussed in the paper, one can look at certain observables which vanish in the kinematic 
limit (like R = ul/ctt in deep-inelastic scattering). Those include <7l, vtt and ctlt discussed 
in Ref. fT2"| . One can also check the slow-Q 2 dependence of certain scaling functions. As a 
first guess, one shall be at least at the deep-inelastic kinematic region, and the recoil nucleon 
shall be at backward angles. It would be interesting to demonstrate experimentally that the 
single-quark scattering mechanism is at work. 

One important theoretical issue is whether a factorization theorem exists for the deeply- 
virtual Compton process. The answer is most likely yes for several reasons. First, the tree 
result we obtained can be easily cast into an operator product expansion. Such an expansion 
is generally believed to be valid independent of external states. Second, the radiative cor- 
rections to the tree process are quite similar to the corrections to 7*7r — > 7, a process known 
to be factorizable at one loop level. Nonetheless, one still has to prove explicitly the fac- 
torization theorem for DVCS, which will be done in a separate publication. A factorization 



proof for a similar process, electroproduction of mesons, has appeared recently |L9|. 

DVCS at small A 2 is especially interesting, because the nucleon form-factor suppression 
is small and because A 2 — > limit is relevant to the spin structure of the nucleon. However, 
due to QED infrared divergences, the Bethe- Heftier process becomes dominant there despite 
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the fact that the DVCS cross section scales like the deep-inelastic scattering cross section. 
Thus one cannot go to too small A 2 . One may get around this to some extent by isolating the 
DVCS and BH interference term through single-spin asymmetry and/or combined lepton- 
antilepton scattering. Theoretical study of A 2 dependence of the OFPDs, in particular, its 
relation to the meson-dominance and the exponential decay law for exclusive cross sections, 
is urgently needed. 

The OFPDs depend on four independent variables: the Bjorken-type of variable £, the 
Feynman-type variable x, the momentum transfer A 2 and the virtual-photon mass Q 2 . It 
is unlikely that one can learn the entire parameter space of the distributions in one kind of 
experiments. As our formulas show, different kinds of experiments are sensitive to different 
combinations of parton distributions and hence they are complementary. [The reason that 
the nucleon helicity-flip distributions contribute to unpolarized scattering is similar to that 
the Pauli-form factor F 2 (A 2 ) appears in unpolarized elastic scattering.] However, to be able 
to study sum rules, one must have accurate data in an extended kinematic region. To achieve 
this, one must have dedicated experiments at a suitable machine like ELFE for extended 
running ||20|| . 

The off-forward parton distributions can also be defined for quark helicity-flip (chiral- 
odd) correlations and for higher twists. DVCS provides one process to access to these 
distributions. There are other processes one can consider to measure them. For instance, 
the diffractive p or J ftp production studied by Brodsky et al. [^TJ can be used to measure the 
off-forward gluon distributions. Recently, Radyushkin has published a paper aimed at this 



direction Jl6|. Thus there is now a new territory to explore the quark and gluon structure 
of the nucleon besides the traditional inclusive (parton distributions) and exclusive (form 
factors) processes. 
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FIGURES 

FIG. 1. Virtual Compton scattering process. 

FIG. 2. QCD diagrams for deeply- virtual Compton scattering: a), the hand-bag diagram, b). 
& c). some 1/Q 2 corrections, d). & e). some radiative corrections. 

FIG. 3. Leading- log evolution of the off-forward parton distributions. 

FIG. 4. Electroproduction of a photon off the nucleon: a), the virtual Compton scattering, b). 
& c). the Bethe-Heitler process. 

FIG. 5. Comparison of the DVCS and BH cross sections at a). CEBAF and b). HERMES 
kinematics. 
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